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We calculate the transmission of electrons and holes between two normal- metal electrodes (N), 
separated over a distance L by an impurity-free superconductor (S) with d-wave symmetry of the 
order parameter. Nodal lines of vanishing excitation gap form ballistic conduction channels for cou- 
pled electron-hole excitations, described by an anisotropic two-dimensional Dirac equation. We find 
that the transmitted electrical and thermal currents, at zero energy, both have the pseudodiffusive 
1 / L scaling characteristic of masslcss Dirac fermions — regardless of the presence of tunnel barriers 
at the NS interfaces. Tunnel barriers reduce the slope of the 1/L scaling in the case of the electrical 
current, while leaving the thermal current unaffected. 
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I. INTRODUCTION 

PseudodifFusive transmission refers to the 1/L scaling 
of the electrical cnirrent transmitted over a distance L 
through a clean sheet of undoped graphene [1]. The same 
1 /L scaling characterizes diffusion in a random potential, 
but now it applies in the absence of any disorder. There 
is a large number of theoretical [2-12] and experimental 
[13-15] studies of this phenomenon, which is understood 
as a general property of masslcss Dirac fermions in the 
limit of vanishing excitation energy. The optical analogue 
in a photonic crystal with a Dirac spectrum has been 
studied as well [16-19]. 

Layered superconductors with a d-wave symmetry of 
the order parameter (notably the high- Tc cuprates [20]) 
form an altogether different system in which massless 
Dirac fermions are known to exist [21-23]. These are 
so-c;alle(l nodal fermions, located in the two-dimensional 
Brillouin zone near the intersections (nodal points) of the 
Fermi surface with lines (nodal lines) of vanishing exci- 
tation gap. Elastic mean free paths I as large as 4 /im 
have been reached in YBa2Cu307_e single-crystals [24], 
much larger than the superconducting coherence length 
^0 — 2nm. It is the purpose of this work to demonstrate 
theoretically the pseudodiffusive \/L scaling of the trans- 
mission through a (i-wave superconductor over the range 
of lengths between and I. This anomalous scahng was 
not noticed in earlier studies of similar systems [25-27] . 

The problem is interesting from a conceptual point of 
view, because it highlights both the differences and simi- 
larities between Dirac fermions produced by a bandstruc- 
ture (as in graphene or photonic crystals) or produced 
by a d-wave order parameter. In undoped graphene, the 
transmitted electrical current I in response to a voltage 
difference V scales as [2, 3] 



h ttL 



(1.1) 



the metal contacts, but small compared to the mean free 
path /. The length L should also be small compared to 
the transverse width W of the graphene sheet (to avoid 
edge effects). Potential barriers (smooth on the scale of 
the lattice constant) at the interfaces between the metal 
contacts and the graphene sheet have no effect on the 
current, because of the phenomenon of Klein tunneling 
[10]. 

For the rf-wave superconductor, we find a transmitted 
electrical current per layer equal to 
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The length L over which the current is transmitted 
should be large compared to the Fermi wave length Af in 



for £,0 <t: L l,W. Here ri_2 € (0,1) are the tun- 
nel probabilities through the potential barriers at the 
two normal-metal-superconductor (NS) interfaces. The 
Dirac equation for nodal fermions is anisotropic [21], with 
different velocities vp and f a parallel and perpendicular 
to the nodal lines. This anisotropy (with vp/va ~ 15 
in YBa2Cu307_e) increases the slope of the 1/L scaling. 
Remarkably enough, the anisotropy does not introduce 
a dependence of the transmitted current on the angle a 
between the direction of the current and the nodal lines. 
The result (1.2) holds generically for any orientation, ex- 
cept for a narrow range of angles of order ^o/L around 
a = (mod 7r/4). 

The tunnel barriers reduce the slope of the 1/L scal- 
ing of the transmitted electrical current (1.2), by a factor 
P1P2/4 for small tunnel probabilities. This does not im- 
ply that the nodal fermions are only weakly transmitted, 
but rather that the transmission probabilities for trans- 
mission as an an electron or as a hole are almost the same 
for Pi 2 <SC 1. Indeed, we find that the electrical shot 
noise power P as well as the transmitted thermal current 
-^thermal (both of which do not depend on the sign of the 
carriers charge) remain finite in the limit ri_2 0. We 
interpret this result in terms of a resonant coupling via 
the nodal lines of the mid-gap states [28, 29] extended 
along the two NS interfaces. We also find, quite surpris- 
ingly, that the thermal conductivity is independent of the 
tunnel probabilities Pi 2- 
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FIG. 1: Geometry to measure the transmission of nodal 
fermions through a d-wave superconductor. A current Ji is 
injected into the superconductor from metal contact Ni (at 
a voltage V) and drained to ground via the superconductor 
(current Is) or via a second metal contact N2 (current I^). 
If the separation L of the metal contacts is large compared 
to the superconducting coherence length ^0, the current I2 is 
predominantly due to transmission parallel to the nodal lines 
a; = or 1/ = of vanishing excitation gap . 



The outline of this paper is as follows. In Sec. II we for- 
mulate the scattering problem and calculate the transfer 
matrix of the nodal Dirac fermions through the rf-wave 
superconductor. The matching of wave functions at the 
interface with the metal electrodes is done in Sec. Ill, 
both for ideal NS interfaces and for interfaces containing 
a tunnel barrier. The transmission matrix of electrons 
and holes follows in Sec. IV. We then apply this result 
to the calculation of transport properties: the electrical 
current (Sec. V), the thermal current (Sec. VI), and the 
electrical shot noise (Sec. VII). We conclude in Sec. VIII 
with a discussion of our results and an outlook. 



II. TRANSFER MATRIX FOR NODAL 
FERMIONS 

A. Anisotropic Dirac equation 

We consider a two-dimensional spin-singlet supercon- 
ductor (S), connecting two normal metal contacts with 
parallel NS interfaces, separated by a distance L. The 
transverse dimension W of the superconducting strip (in 
the x-y plane) is assumed to be large compared to L, in 
order to avoid edge effects. The order parameter A(fc) 
is assumed to have d^y symmetry: it vanishes for wave 
vectors along two nodal lines, which are taken to be the 
X and the y axis. All our results also apply to dx^-y^- 
superconductors, for our purposes, a simple 7r/4 rotation 
relates the two systems. To be specific, the x-y plane can 
represent a single Cu02 layer of a cuprate superconduc- 
tor [20], with the [100] direction at an angle 7r/4. 

Low-energy excitations in the superconductor are 
found in the Brillouin zone near the four intersections 
(±fci?,0), (Ojzkkp) of the Fermi surface with the nodal 
lines of the order parameter. (These nodal points are la- 
beled A,B,C,D in Fig. 2.) Around these points, both 




FIG. 2: Ellipsoidal equal-energy contours of low-energy ex- 
citations in the Brillouin zone of a superconductor with dxy 
symmetry. Long and short axes have ratio vf/va- The con- 
tours are centered at the four nodal points (solid dots), where 
the order parameter vanishes on the Fermi surface. The nor- 
mal n to the NS interfaces is indicated. The dashed line, 
displaced from the nearest nodal point by q, indicates points 
of constant wave vector component parallel to the interface. 

the pair potential A(fe) and the kinetic energy can be lin- 
earized: the dynamics of the nodal fermions is governed 
by an anisotropic Dirac equation [21-23]. For example, 
near node A at (kp,0) this can be written in the form 

or more compactly with the help of Pauli matrices, 

-i[vFaA + vi:^a:,dy\^ = e-^. (2.2) 

We have set h to unity, restoring units in the final ex- 
pressions. The spinor — (^'e,^';i) contains the en- 
velope wave functions of electron and hole excitations 
(slowly varying on the scale of the Fermi wavelength 
\p = 2Tr/kp). The Fermi velocity vp is larger than the 
velocity va ~ Aq /hkp by a factor of order £,o/Xp (with 
^0 = fivp/IS.Q the superconducting coherence length), 
which is in the range 10-20 for cuprate superconductors. 
The equal-energy contours in the Brillouin zone of the 
nodal fermions thus have an elongated ellipsoidal shape, 

e{5k) = ^J{vp5kxY + {vA5ky)\ (2.3) 

as a function of the displacement 5k of the wave vector 
from the nodal point. 

B. Transfer matrix 

Since the system is translation invariant along the NS 
interfaces, the component of the wave vector along these 
interfaces, q = —Sk^ sin a+Sky cos a, is a conserved quan- 
tity. Here a is the angle between the normal to the NS 
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interface and the nodal line pointing to node A, which 
we restrict to— 7r/4<a<7r/4 without loss of generality. 
Moreover, since mirror reflection along the NS interface, 
followed by the transformation A(fe) — ^ — A(A;), while 
leaving all the other parameters unchanged, maps a on 
—a, we can further restrict a to < a < 7r/4. In all our 
formulas, to obtain the corresponding formulas for —a, 
replace q by —q and f a by —va- 

We write "^{r) = \l'(s)e*«'*', with s e (0,i) the co- 
ordinate perpendicular to the NS interfaces and s' the 
coordinate parallel to them. We substitute ^(r) into 
Eq. (2.1) and find that the spinor ^'(s) satisfies the wave 
equation 



+ Ji(z]*(s) = £*(s), 



(2.4) 



where dg is differentiation perpendicular to the NS inter- 
face, and J and Ji are the operators of particle current 
perpendicular and parallel to the NS interface, 

J = vpCTz cos a + vacTx sin a , (2-5) 

Ji = vao'x cos a — vpCTz sin a. (2.6) 

We note that the operator J squares to a scalar, its mag- 
nitude giving the particle velocity Va perpendicular to 

the NS interface: 



v^ = = Vp cos^ a + v"^ sin^ a. 



(2.7) 

To solve Eq. (2.4), we multiply it by J/v"^ and rear- 
range to obtain 



(2.8) 



with 



^0 = 9 



sin 2a vi 



. VpVA 



J-iq^^ay. (2.9) 



2 vl - 

The solution to Eq. (2.4) can then be written as 

*(so + s) = A^s*(so); Ms = exp[Uos], (2.10) 

where the second equation defines the transfer ma- 
trix Ms- As expected, the particle current J is con- 
served by Eq. (2.4): [*^(s) J«'(s)] = *^(s)[-M^J-F 
iJAo]'^{s) = 0. 



III. WAVE MATCHING AT THE NS 
INTERFACES 

At the two NS interfaces the coupled electron-hole ex- 
citations in the superconductor are converted into uncou- 
pled electrons and holes in the normal metal. We thus 
need to match, at s = and s = L, the envelope wave 
functions ^' = (5',., 5';j) of the nodal fermions in S to the 
Bloch wave functions $ = {'S>,.,^h) of free fermions in 
N. This is similar to the matching of Dirac equation to 
Helmholtz equation considered in the context of trans- 
mission through a photonic crystal [16]. Translational 
invariance parallel to the NS interfaces requires that the 
coupling conserve the wave vector component q paral- 
lel to the interfaces. Particle flux conservation imposes 
further constraints, as we determine here. 



A. Particle flux conservation at the NS interface 

At the surface of the superconductor, the order param- 
eter A attains its bulk value over a short length scale, the 
healing length Iq. The two-component wave function on 
the S side of the interface (at s = Iq) can be linked to 
that on the N side (s = 0) by an interface matrix A4ns, 
defined by 



*(W = ^ivs*(0). 



(3.1) 



In the normal metal, the operator of particle flux per- 
pendicular to the NS interface can be written as 



JN = VnCTz, 



(3.2) 



with vn possibly different from vpcosa because of a 
Fermi energy mismatch. The requirement of particle flux 
conservation reads 

*(0)ti;jva,*(0) = *(0)tM]v5JMjvsf (0). (3.3) 

To derive the most general form of the interface matrix 
fulfilling this requirement, notice that a unitary rotation 
through angle 6, where 



9 = arctan 



VA 



tana 



transforms J into up to a scalar factor: 

Vaf^z = exp[i6'CTj,/2] Jexp[— i6'(Tj,/2]. 
This allows us to write the interface matrix as 



(3.4) 



(3.5) 



Mns = .l^exp[-ieay/2]Mo, (3.6) 

V 

where A4q is a 2 x 2 matrix fulfilling a generalized uni- 

tarity condition. 



(3.7) 



Eq. (3.7) restricts Mq to a three-parameter form 

Mo = e'^''''''e^^''^e'^''''' (3.8) 

(ignoring an irrelevant scalar phase factor), with arbi- 
trary real parameters Px,Py,Pz- To understand better 
where the nontrivial interface matrix arises from, and to 
show that we may set (3x = Py — Pz = 0, we have to 
extend the Dirac equation (2.4) to the interface layer, 
where va varies in space. This is done in Appendix A. 

So far we have considered only intranode scattering at 
the NS interface. We refer to such an interface as an 
"ideal interface" . A nonideal interface contains a tunnel 
barrier, which introduces internode scattering. We will 
consider the transfer matrices through the rf-wave super- 
conductor for both cases in the next two subsections. 
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B. Transfer matrix with ideal NS interfaces 

The complete transfer matrix for a strip of rf-wave su- 
perconductor with ideal NS interfaces reads 



Xideal = Mj^^MlMns, 



(3.9) 



where Alz, is the Ms from Eq. (2.10) with s ^ L, describ- 
ing propagation inside the superconductor, and Mns 
from Eq. (3.6), with Mo = 1, describes an NS interface. 
Upon substitution, we obtain 



■Mideal = e 



+ 



cosh.[Ka{q)L\ 
sinh[KQ,(g)L] 



with the definitions 



ba{q) = qL- 



sin 2a 



(3.10) 

(3.11) 
(3.12) 



Notice, how as a result of accounting for the two 
NS interfaces ~ the transfer matrix has simphfied from 
that of Eq. (2.10). The change is that the particle flux 
operator J in Eq. (2.9) is replaced by Vacr^ in Eq. (3.10). 
Also note that the determinant of the transfer matrix 
has norm one, DctA^idoai = e^"'^"'^^ as required by the 
generalized unitarity relation 



M- 



(3.13) 



which holds for any transfer matrix as a consequence of 

particle current conservation. 

To appreciate the effects of the Dirac cone anisotropy, 
we can perform a linear transformation on our system 
to obtain one with an isotropic Dirac cone: contraction 
along the nodal line by a factor v^/va, and expansion 
perpendicular to it by a factor vp/va- The dispersion 
of the new, isotropic Dirac cone has a single velocity pa- 
rameter V = vpVA/va- The superconducting strip is de- 
formed by the transformation: its width W is unchanged, 
but its length L becomes 



L 



VpVA 



V 



(3.14) 



an effective propagation length we define here for later 
use. 

The matrix (3.10) derived above is the transfer ma- 
trix for nodal fermions near point Ua = (A;ir,0,0) on the 
Fermi surface, with q = qA = {k — fe^) • (z x n) the trans- 
verse wave vector component relative to fc^. Similarly, 
the transfer matrices near each of the four nodal points 



can be written as 
Ma = e^^«(-) expf ^a. + ^-^^a^) , (3.15a) 

= e-^w.-o(..)expf -^a. - ^^^a.V 

(3.15b) 

Mc = e^^»(-)expf-^a. + ^^^^.,1 (3.15c) 
Md = e"'*'''/^-°(«")exp <jz - ^\ 



7r/2-a 



(3.15d) 

The basis at each nodal point is the same spinor (<i>e, <i>/i), 
but the electron states $e are "right-movers" (propagat- 
ing from Ni to N2) at nodal points A or B and "left- 
movers" (from N2 to Ni) at nodal points C and D. 



C. Nonideal interfaces 

The complete Fermi surface of the normal metal (N) 
might differ in many ways from that of the superconduc- 
tor. However, when we study transport near a specific 
nodal point, due to transverse momentum conservation, 
we can effectively reduce the Fermi surface to the two k 
points where transverse momentum has the same value 
as at the nodal point. These two k points in N each cou- 
ple to different nodal points in S, for example to nodal 
points A and C in Fig. 2. A nonideal NS interface cou- 
ples different nodal points, by reversing the component 
of the momentum perpendicular to the interface. Such 
internode scattering may be caused by an insulating layer 
at the NS interface, or it may result from the Fermi veloc- 
ity mismatch between N and S. Note that only internode 
scattering is possible in the absence of superconducting 
order — any intranode scattering has to happen inside 
the superconductor. 

We will generically describe a nonideal NS interface 
by a tunnel barrier, with tunnel probability F (which we 
take mode independent for simplicity). For \a\ < £,o/L, 
the tunnel barrier couples electrons neax nodal points A 
and C. The transfer matrix Mr{so) for a tunnel barrier 
at position sq, defined by 



e,A 
e,C 



•.,C 



(3.16) 



has the form 

, . , , [if 1 e-''^('*°)yr^^ 

Mt{so) = y ^ [e'Hso)^iZIT 1 

(3.17) 

with (j){s) — 2kFS cos a. 

The tunnel barrier at Sq also couples holes near nodal 
points A and C, with transfer matrix 



^'-^^ = A^* (so) ^^'^.c: 



h,A 



(3.18) 
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(The basis states are chosen such that the upper com- 
ponent is a right-mover and the lower component a left- 
mover.) 

Finally, we can write down the full transfer ma- 
trix of the superconducting strip, in the basis 
{^e^A,^h,A,^e,c,^h,c), Including nonideal contacts 
with tunneling probabilities Fi at s = and F2 at s = L. 
It is obtained by matrix multiplication. 



M = V 











V 



M*r,iL) 



Mc 











(3.19) 



with V a unitary matrix that switches bases from 

{^e,A,^e,C,^h,C,'^h,A) {'^>e.A,^h,A,'^e,C,'^h,c)- 



V 



/I ON^ 

1 

10 

Vo 1 0/ 



(3.20) 



If 7r/4 — a < £,q/L, the nodal point A is coupled to the 
nodal point D, so the above formulas still hold, with C 
replaced by D. 

If both a £,o/L and 7r/4 — a ^ ^o/L, the tunnel 
barriers at the interfaces do not couple nodal point A 
to any other nodal points. Since we assume ^o/L ^ 1, 
this case of misaligned nodes is the generic case. In that 
case, A^c is to be replaced by the singular transfer ma- 
trix A^Andrccv Corresponding to Andreev reflection with 
reflection amplitude — «, 



■M Andreev = Am 6" 



(3.21) 



Since Andreev is also the qc — > — oo limit of Mc{qc) in 
Eq. (3.15c) (up to an irrelevant phase factor), Eq. (3.19) 
is valid as it stands for misaligned nodes as well. 



IV. TRANSMISSION AMPLITUDES 

A. Ideal interfaces 

Referring to the geometry of Fig. 1, a scattering state 
(for a given value of q) has the form <i>(0) — (l, r^g) at the 
normal side of the left NS interface and <&(i) — (tee,0) 
at the normal side of the right NS interface. The com- 
plex number r^e is the amplitude for Andreev reflection 
(from electron to hole) and the complex number tee is 



the amplitude for electron transmission. We calculate 
this transmission amplitude using the relation 

tee=([A4^']ll)"'=([M]ll)"\ (4-1) 

where the first equality follows from $(0) = A4^^^{L), 
and the second equality from particle current conserva- 
tion, Eq. (3.13). 

Substitution of Eq. (3.15a) gives the expression 



pi<t>a(qA) 



cos\i[Ko,{qA)L] 



iesmh[Ka{qA)L\ 

VaKaiqA) 



(4.2) 



Nonideal interfaces 



For nonideal interfaces we have to consider both the 
transmission amplitude tge from electron to electron and 
the transmission amplitude the from electron to hole. It 
is convenient to define the 2x2 transmission matrix 



t = 



the thh 



(4.3) 



which contains also the transmission amplitudes teh and 
thh from hole to electron and from hole to hole. This 
matrix t is a 2 x 2 subblock of the 4x4 unitary scattering 
matrix S, which we derive in Appendix B. 

To obtain t from the 4x4 transfer matrix M , we make 
a change of basis from the basis {^e,A, ^h,A, ^e,C, ^h,c) 
used in Eq. (3.19) to a basis {^e,A,^h,c,^e,c,^h,A) in 
which the upper two components are right-movers and 
the lower two components are left-movers. The change 
of basis is carried out by the unitary matrix 



W 



/I 0^ 

1 

10 

\o 1 oy 



(4.4) 



We can then follow the same reasoning as in the previous 
subsection, to conclude that t is determined by the 2x2 
upper-left block Xn of WMW^^, 



xlt 



1, 



(4.5) 



cf. Eq. (4.1). 

Substitution of A4 from Eq. 

algebra. 



(3.19) gives, after some 
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. ^ jViT^y/^ / {Mc)22 + eM^MA)22VT^iVT^ -{Ma)i2VI^i - e-"*'<-'^HMc)i2VT^2\ 
Z \-{Mc)2iVl - Ti - e^'^(^H^A)2ivT^ {Ma)ii + e-'*^^'>{Mc)iiVT^iVT^ J ' 

(4.6) 

Z = ^/l-Tl^/l-T2 (e-'^^^^DetMc + e^'I'^^^'DetMA) + {Ma)ii{Mc)22 

+ iMA)22iMc)llil - ri)(l - T2) - (Ma)i2(Mc)2i(1 " Ti) - {MA)2liMc)l2il - T2). (4.7) 

I 



V. ELECTRICAL CURRENT 
A. Ideal interfaces 

Turning now to observable quantities, we will work in 
the linear response regime V ^ 0, when the transmission 
amplitudes may be evaluated at the Fermi level (e = 0). 

The current (per layer) transmitted into metal con- 
tact N2 through nodal point A is obtained by integrating 
the transmission probability |ieeP over qa, 



GoV 



W 
2n 



dQA \te 



(5.1) 



(The conductance quantum Go = 2e^//i includes a 
twofold spin degeneracy.) The integrand decays expo- 
nentially for > v'^/vfVaL 0^ {S,o/L)kF- For i > ^0 
the effective integration range is much smaller than kp 
and may be extended to ±00. Substituting Eq. (4.2) (for 
e = 0) we arrive at 



GoV 



W 



(5.2) 



L ttvfVa 

As expected, the conductance of a single nodal point has 
the same form as that of a single valley in a graphene 
strip, with L replaced by the effective propagation length 
L„ of Eq. (3.14). 

The current transmitted through nodal point B is 
given by the same formula with Va replaced by Ua-7r/2- 
Because of the identity 



2,2 2,2 



(5.3) 



the total current I2 = I2 + I2 becomes independent of 
a. The conductivity aideai = {l2/V){L/W) per layer for 
the case of ideal NS interfaces is then equal to 



Cideal 



(5.4) 



mipVA 

As discussed in Sec. VIII A, Eq. (5.4) differs [by a factor 
1 + {va/vf)'^] from the bulk electrical conductivity of 
Refs. [21, 22]. 



B. Nonideal interfaces 

For nonideal NS interfaces, tunnel barriers couple the 
nodal points, and the calculation of the current 1 2 be- 
comes more involved. In this Section we treat the generic 



case of misaligned nodal points. The case of (perfectly) 
aligned nodal points is considered in Appendix C. 

We first calculate the current through nodal point A. 
As discussed in Sect. IIIC, we can substitute Mc with 
A^Andreev of Eq. (3.21), and using Eq. (4.6) we obtain 
the transmission matrix (at e = 0) 



ei't'o.MZA 
1 



where the denominator Za has the form 



(5.5) 



Za = Vri(2 - ri)r2(2 - T2) cos\i[L^{qA - ^peak)]- 

(5.6) 

Here La is the effective propagation length (3.14), while 
9peak is the transverse wave number defined by 



5peak 



2L, 



In 



2 - Fi 2 - Fs 



(5.7) 



Both fee and th^ are peaked at qa = 9peak- This peak 
momentum lies at the nodal point (^pcak = 0) only for 
ideal interfaces. In the presence of tunnel barriers the 
sign of gpeak is such that the order parameter has opposite 
sign at the two intersections of the line qa = 9peak with 
the Fermi surface. 

Integration over qa of electron current minus hole cur- 
rent gives the net (electrical) current. 



-'2 



GoV 
GoV 



W 
W 



dQA [\{tA)e 



\{tA) 



he\ 



L tvvfVa 2 - Fi 2 - F2 



(5.8) 



Similarly, for the current through nodal point G we 
take the limit — > 00 of Eq. (4.6) and obtain the trans- 
mission matrix 



V^2 



e^MQc) Zc 



2 I 



(5.9) 



Zc = Vri(2 - ri)r2(2 - T2) COSh[L„(gA + apeak)], 

(5.10) 
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and then the current 



-'2 



-GoV 



L Tivpv^ 2 — Fi 2 — r2 



(5.11) 



Note the minus sign in the formula for I2 . The current 
has opposite sign to that at nodal point A, since here 
holes rather than electrons tunnel across the system to 
contact 2. 

The total current (per layer) through nodal points A 
and C becomes 



G.V 



W 



L ttvfva 2 — Fi 2 — r2 



(5.12) 



Comparison with Eq. (5.2) reveals that each tunnel bar- 
rier changes the sum of the current transmitted through 
a nodal point and the one opposite to it in momentum 
space, its time-reversed partner, by a factor of r/(2 — F). 
As in the case of ideal NS contacts, the pair of nodal 
points B and D contribute a same amount, but with Va 
replaced by v^j^-a. The a-dependence again drops out 
of the total current h = + 1^ + 1^ + 1^ ■ For the con- 
ductivity u = {l2/V){L/W) per layer we finally obtain 



2-ri2-r2 



(5.13) 



VI. THERMAL CURRENT 

The conductivity (5.13) vanishes in the weak tunneling 
limit Fi, F2 0, because the electron and hole contribu- 
tions to the electrical current I2 then become equal but 
of opposite sign. Electrons and holes contribute with the 
same sign to the thermal current, 

W 

/thermal = ioGoTJr— y dq [\tee\^ + \the\^] , (6.1) 

with Lq = ■K^k'^/Ze^ the Lorenz number. The thermal 
current flows from contact Ni at temperature T + 5T 
into contact N2 at temperature T. (Eq. (6.1) requires 
5T <C T and T sufficiently small that the transmission 
amplitudes may be evaluated at the Fermi energy e = 0.) 

We consider the (generic) case of misaligned nodes. 
Substitution of the expressions for t from Sec. VB, and 
summing over the pair of nodal points A and C, we find 
that 



tA , jC 

-'thermal ^ 



thermal 



L0G0T6T 



W vl 



(6.2) 



Quite surprisingly, this turns out to be independent of 

the tunnel probabilities Fi and F2. The total thermal 
current (per layer) also includes contributions from the 
nodal points B and D, and is just as the electrical 
conductivity - independent of the angle a: 



/thermal = LqGoTST- ^ 



(6.3) 



As discussed in Sec. VIII A, the thermal conductivity k = 
{Ithermai/ ST) (L/W) extracted from Eq. (6.3) coincides 
with the bulk thermal conductivity of Ref. [22]. 



VII. SHOT NOISE 

The zero- frequency noise power of time dependent elec- 
trical current fluctuations 5/2 (i) measured in contact 
number 2, 



/oo 
dtSl2iO)Sh{t), 
-00 



(7.1) 



is given in terms of the transmission matrix elements by 
the general expression [30] 

P22 = GoeV^ J dq [\t,,f{l - |iee|') 

+ \thef{l - \thef) + 2\the\''\tee\''] ■ (7.2) 

As with the conductance, we work in the linear response 
regime, so the transmission matrix is to be evaluated at 
e = 0. 

We restrict ourselves to the case of misaUgned nodes 
and substitute the expressions for t from Sec. VB. The 
integral over q contains four separate contributions, from 
q near nodes A, B, C, and D. The total result (per layer) 
is 



GoeV 



L itvfva 

i2(2-ri)2(i 



F2) + 8(i-Fi)f: 



3(2-ri)2(2-F2)2 

The Fano Factor F = P22/e/2 is given by 

12(2 - Fi)2(l - F2) + 8(1 - Fi)F2 + F^r 



F = 



3rir2(2-ri)(2-r2) 



(7.3) 



(7.4) 



In the ideal limit Fi, r2 1 we find a Fano factor F = 
1/3, three times smaller than the value F = 1 associated 
with a Poisson process. As discussed in the context of 
graphene [3, 14, 15], this is the same one-third reduction 
as in a diffusive metallic conductor and is a hallmark of 
pseudodiffusive transmission. 

In the weak tunneling limit Fi, r2 — > the noise power 
remains finite, 



lim P22 = GoeV ^ 
ri,r2^o L 



Wv% + v% 



(7.5) 



while the electrical current vanishes, I2 oc F1F2 0. The 
electrical current fluctuations therefore become large rel- 
ative to the time-averaged current in the presence of tun- 
nel barriers. This is discussed in the context of resonant 
tunneling through midgap states in Sec. VIII B. 
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VIII. DISCUSSION 



resulting from that work, 



Comparison with bulk electrical and thermal 
conductivities 



The electrical current I2 and thermal current /thermal 
that we have calculated describe transniission of electrons 
and holes over a finite length L of a clean d-wave su- 
perconductor. Earlier work [21, 22] calculated the elec- 
trical and thermal conductivities cto and kq of a dis- 
ordered infinite system. These are in principle differ- 
ent systems, but we can still compare them by formally 
converting the currents through the finite system into 
bulk conductivities by means of ao = {l2/V){L/W) and 

no = (/thermal/5T)(L/W^). 

The thermal conductivity obtained in this way from 
the finite-system thermal current (6.3), 



Ko = LqGqT 



TTVfVa 



(8.1) 



is the same as the bulk thermal conductivity of Durst 
and Lee [22]. The results for the electrical conductivity 
differ, however. The bulk result [21, 22] 



ao = Go 



VF 



(8.2) 



differs from the finite-system result (5.4) — even if we 
assume ideal NS interfaces. The difference between the 
factors vf/va in Eq. (8.2) and (vp + v'^)/vfva in Eq. 
(5.4) is small in practice (because vf ^ va), but the 
difference does illustrate that these are different systems. 



B. Interpretation in terms of resonant tunneling 
through midgap states 



\tehf = \the\'^ = \thhf 



(8.3) 

with the results from Sec. VB in the tunneling limit 
Ti,T2 « 1, 



1 



4 cosh [ia(a - 9peak)] 



We have defined 



This comparison leads to the identification 



2e^«« 



(8.4) 



(8.5) 



(8.6) 



Resonant tunneling, with all transmission probabili- 
ties equal to 1/4, occurs when q = (/peak, hence when 
-E'edge = VSEqSEl (tunnel splitting of the midgap states 
equal to tunnel broadening). Because transmission from 
electron to electron and from electron to hole happens 
with the same probability (to leading order in ri.r2); 
the transmitted electrical current vanishes in the limit 
of small r. The thermal current /thermal and electrical 
noise P22 remain finite, because [teej^ and Ithe]"^ con- 
tribute with the same sign to these quantities. 

This interpretation explains the finite small-F limit for 
P22 and /thermal, but it does not explain why the thermal 
current (6.3) turns out to be completely independent on 
the values of Fi and That remains a surprising re- 
sult of our calculation, for which we have no qualitative 
explanation. 



Outlook 



We have found that tunnel barriers at the NS inter- 
faces reduce the transmitted electrical current, but not 
th{^ thc^rmal current nor the electrical noise. This result 
has a natural interpretation in terms of the midgap states 
at the NS interfaces. Midgap states are zero-energy edge 
states of the rf-wave superconductor, which exist at mo- 
mentum q along the edge if the order parameter has op- 
posite sign at the two intersections of the line of constant 
q with the Fermi surface [28, 29]. The midgap states at 
the two NS interfaces have a small overlap, and there- 
fore acquire a nonzero energy ±i?edgo (tunnel splitting). 
Moreover, the coupling to the metal electrodes at s = 0, L 
introduces partial widths 6E0, SEl of the midgap states 
(tunnel broadening). 

Tunneling through a pair of midgap states was stud- 
ied in Ref. [31], in the context of Majorana bound 

states (which are a special type of nondegenerate midgap 
states). We can compare the transmission probabilities 



We have shown how ballistic transport through a 
clean rf-wave superconductor (such as single-crystal 
YBa2Cu307_e) has features in common with graphene 
[1]: a pseudodiffusive 1/L scaling of the electrical current 
transmitted over a distance L, and a 1/3 suppresssion of 
the electrical shot noise with respect to the Poisson value 
of uncorrelated current pulses. These effects have been 
observed in graphene [13 15] and it would be interest- 
ing to search for them in the high-Tg cuprates. The 1/L 
scaling should persist, with a modified slope, in the pres- 
ence of tunnel barriers at the NS interfaces, and in the 
case of the thermal current we find that even the slope is 
independent of the tunnel barrier height. 

There are more areas of correspondence between mass- 
less Dirac fermions in d-wave superconductors and in 
graphene, in addition to the pseudodiffusive transport 
studied in this work. We mention two such effects, as 
directions for future research. 
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• In graphcnc an electrostatic potential can displace 
the Fermi level away from the Dirac point of van- 
ishing density of states. In the d-wacve supercon- 
ductor the supercurrent velocity Vs enters into the 
Dirac equation (2.2) as a scalar term oc ctq [32], 
and therefore has the same effect of displacing the 
Dirac point relative to the Fermi level. There is one 
curious difference with respect to graphene: the d- 
wave superconductor has two pairs of valleys and 
the Dirac point can be displaced independently in 
each pair (relative to the same Fermi level). With 
reference to Fig. 2, the component of Vs in the x- 
direction acts on valleys at the nodal points A and 
C, while the component in the ?/-direction acts on 
those at B and D. 

• While the role of an electrostatic potential in 
graphene is played by the supercurrent, an electric 
field in the d-wave supcrconciuctor plays the role of 
a magnetic field in graphene. If a sufiiciently strong 
electric field could be induced in a thin-film cuprate 
superconductor, it might be possible to see effects 
analogous to the effects of Landau level quantiza- 
tion in graphene [33]. 
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s-dependent and gets a new term: 



A{s) = Aois) 

v'^{s) sin a 



{ayVp cosa — ivA{s)sma) . (A3) 



Since qloVFVA/Va — l/kp ^ 1 (in the relevant range 
of g's near the nodal point), the integral of Aq over the 
interface layer is <C 1 and may be neglected. Then A{si) 
commutes with A{s2) for < Si,S2 < ^O) and therefore 
we can simply integrate Eq. (2.8) over the interface layer: 



Mns = exp 



= exp 



pit) 

i / A{s)di 
Jo 



{ay — iu) du 



vp cos a 



exp[—i6ay/2], 



(A4) 



with 6 = arctan [(iia/^f) tana] as defined in Eq. (3.4). 
The result agrees with Eq. (3.6) with f]^ = Py = f3^ = 0, 
and Vn = vpcosa. The Fermi velocity mismatch con- 
tributes an additional factor \/vnJvf to the interface 
matrix, and in addition may cause internode scattering 
(as detailed in Sec. IIIG). 



APPENDIX B: FULL SCATTERING MATRIX 



APPENDIX A: NS INTERFACE MATRIX 

In Sec. IIIB we derived the most general form of the 
transfer matrix of an NS interface, consistent with the 
requirement of particle flux conservation. The result in 
Eq. (3.6) has three undetermined parameters j3,j,, j3y, and 
[3z- Here we calculate the interface matrix by solving the 
Dirac equation in the interface layer and determine these 
unknown parameters. The interface layer is the region 
where the order parameter increases from to its bulk 
value, over a healing length Zo (wliic;li is typic;ally of the 
same order of magnitude as the coherence length ^o)- 

As discussed in Ref . [34] , in order to preserve Hermitic- 
ity, the Dirac equation (2.2) needs to be supplemented by 
terms containing the spatial derivatives of v/^: 

-i [vpCTzdx + VAC^xdy + {dyVA)'Jx/'A * = (Al) 

We assumed that the phase of A is constant, and set it to 
(without loss of generality), thus ?;a is real throughout, 

and is only a function of the distance ,s from the NS 
interface. An eigenstate ^{s) of momentum q parallel to 
the NS interface satisfies 

[-Uds + Jiq - Wa sin(a)£7^/2] '^{s) = s^{s), (A2) 
with the derivative of f a denoted by the shorthand v'^ = 
dsVA- Accordingly, the matrix A in Eq. (2.8) becomes 



In Sees. IV and V we have calculated the 2x2 trans- 
mission matrix t, which is the quantity we need for the 

transport properties considered. For reference, we give 
here the full 4x4 scattering matrix, 

^ =(;:;), (Bi) 

containing the 2x2 transmission matrices t (from left to 
right) and t' (from right to left), as well as the reflection 
matrices r (from left to left) and r' (from right to right). 
These matrices can be obtained from transfer matrix A4 
by constructing the four 2x2 sub-blocks Xij, 

and then evaluating 

r = -X^2^2\, r' = Xi2X^2^, 

tt=XiY, t' = X^^\ (B3) 

cf. Eqs. (4.4) and (4.5). 

We restrict ourselves to e = and misaligned nodes. 
Near node A we flnd the reflection matrices 
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, ^ J_ fe-'t'(i^^e^-'i^V^^^{2 - Ti) - iT^ie^-'^^ - Ticosh L„«a)\ 



(B4) 
(B5) 



The transmission matrix Ia is given by Eq. (5.5) and = ayt\ay. The resulting scattering matrix (Bl) is unitary, 
SS^ = 1, as it should be. 

Similarly, near node C we find tc given by Eq. (5.9), t'^ = ayt^ay, and the reflection matrices 



rc 



_e-i.9cyr^rT(2-r2) 

1 /e-*'^(^)e-^-«^v'n^(2-ri) -ir2{e-^-i'^ +rismhLo,qc) 



-e-^°9^^-yr^TT(2-r2 

-iri(e"-^°«^' +r2sinhL, 



<lc) 



^ Zc\ -ir2(e-^°9c - Ti coshL^qc) e'^^^h-^""^ ^T^^i2 - Ti) 

r 



(B6) 
(B7) 



APPENDIX C: CONDUCTANCE FOR ALIGNED 
NODAL POINTS 

1. Alignment of nodes A — C. 



and D remain misaligned. Restricting ourselves again to 
£ = 0, we may put qa = qc = q, 4>a = 0, Ma = Mc, 
and La = Lva/vf = Lq in Eq. (4.6). The result is 



For I a I <^ ^,o/L the two nodal points A and C line 
up with the normal to the NS interface, while nodes B 



. _ {TxTiY''^ / [1 + e^^'^^-^yr^^v^n^] cosh Log i[v^r^^ + e-2''=^F-^yr^^] sinh L^q 



Z^C V-«[^r^r7+e2*'=^^\/r^T^] sinh Log [1 + e-^^'^^^Vl^TTyr^I^] coshLog^ ' ^ ' 

Zac = 2\/r^v^r^cos(2fcj.L) + 2 - Ti - Fa + TiFa cosh^ Log. (C2) 

The current 12*^ through the aligned nodes A and C follows from 

lf' = GoV-jjq{\t,,\^-\tU), (C3) 

1^ |2 _ U2 I ^^1^2 

I eel \he\ {2 -Ti){2 -T2) + ^s/l^is/l^2COs{2kFL) +V^V2COsh{2Loqy ^ ^ 



For the total current I2 we add the contribution from the tunnel barriers, 
(strongly) misaligned nodes B and D, 



rmin _ ^ artanhx , ^^a\ 

r.^^^ + Gol/^^^i ^. (C5) ' -^°^L4^^U^~^^J' ^^'^ 

L 7rt;F2-Fi 2-F2 W 1 f vp va \ 



As shown in Fig. 3, the current I2 oscillates as a func- 
tion of fc^L, between minima /™'" at fc^L = (modTr) 
and maxima 7™**^ at fcirL = 7r/2 (modTr). (Similar oscil- 
lations were found in Ref. [25].) Simple expressions for 
these two values follow for the case Fi = F2 = F of equal 



with abbreviation x = 2(2 - F)-Vl - F. For F = 1 we 
recover the ideal limit I^'" = I^^ = aid^^iVW/L. For 
F <C 1 we have instead LJ^^^ = GQV{W/TrL){vF/vA), 
/Jii" = /fia^ X iF2|lnF|. 
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2. Alignment of nodes A — D and B — C 



FIG. 3: Dependence on the separation L of the NS interfaces 
of the current I2 into contact N2, for the interface orientation 

Of = of aligned nodes A and C. Calculated from Eqs. (C3)- 
(C5) for parameters Fi = r2 = 0.3, vf/va = 10. 




For |q!— 7r/4| <C £,i)/L- nodal points A—D and B — C arc 
pairwise aligned with the normal to the NS interface (see 
Fig. 4). The transmission matrix Iad through nodes A — 
D is given by Eqs. (4.6) and (4.7) with Mc replaced by 
M.D- Similarly, for the transmission matrix tsc through 
nodes B — C we should replace AdA^y -Mb- 



FIG. 4: Same as Fig. 2, but now for an angle a = tt/4 between 
the normal n to the NS interface and the lines x = 0, y = 

of vanishing order parameter. For this orientation the nodal 
points A — D and B — C are pairwise aligned with n (dashed 
lines), so that they are pairwise coupled by a tunnel barrier 
at the interfaces. 



Considering first the transmission through nodes A — 
D, we see from Eq. (3.15) that Md = ■M'^^ at \a\ = 7r/4, 
qA = qo = Q- Restricting ourselves to £ = 0, we find 



^ 2(rir2)i/^ / (1 + e"/'('')yr^yr^) coshL„<z -ze-^'^^^^(yr^- e*'^(«)^/^^) sinh L^q\ 
AD gi0,/4(9A)^^^ i^-i(e*'^(«Vl - r2 - - Ti) sinhLag e-*'=^-^^(e*'^(«) + yT^^yT^^) coshL„g /' 

(C8) 

^AB = rir2 + 4v/r^V'r^cosV(g) + (2-ri)(2-r2)cosh2L„g, (C9) 

with La = 2LvfVa/{v% + v\) and the auxiliary function 'tjj{q) = 2Lq{vp — v'^)/{vp + v^) + kpL^/^. 
The current through the aligned nodes A and D follows from 

= GoV- J Jq{\t,e\' - \tl\), (CIO) 
„ ,2 „2 , OF r (2 - ri)(2 - r2) + [T,T2 + 4 VT^iVT^ cos j^jg)] cosh(2L^g) 

[rir2 + 4vr^^vr^T^ COS V'(g) + (2 - ri)(2 - r2) cosh(2L«5)] 

I 



The contribution from the aligned nodes B and C is iden- tical, so the total current becomes I2 = 2/3 



AD 
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For ideal interfaces (Fi = F2 = 1), we recover the re- 
sult I2 = aidea.\VW/L. In the presence of tunnel barriers, 
I2 again oscillates as a function of L, see Fig. 5. 



7r/2 TT 

kFL/V2 



FIG. 5: Same as Fig. 3, but now for an interface orientation 
a = ir/4 of pairwise aligned nodes A—D and B—C, calculated 
from Eq. (CIO) for parameters Fi = r2 = 0.3, vf/va = 10. 
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